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D ■ ABSTRACT 

Cn| : We present an F606W-F814W color -magnitude diagram for the Draco dwarf 

spheroidal galaxy based on Hubble Space Telescope WFPC2 images. The lumi- 
nosity function is well-sampled to ~ 3 magnitudes below the turn-off. We see 
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lO . no evidence for multiple turnoffs and conclude that, at least over the field of the 

! view of the WFPC2, star formation was primarily single-epoch. If the observed 

O 
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number of blue stragglers is due to extended star formation, then roughly 6% 
(upper limit) of the stars could be half as old as the bulk of the galaxy. The color 
difference between the red giant branch and the turnoff is consistent with an old 
population and is very similar to that observed in the old, metal-poor Galactic 
globular clusters M68 and M92. Despite its red horizontal branch, Draco appears 
to be older than M68 and M92 by 1.6 ± 2.5 Gyrs, lending support to the argu- 
ment that the "second parameter" which governs horizontal branch morphology 
must be something other than age. Draco's observed luminosity function is very 
similar to that of M68, and the derived initial mass function is consistent with 
that of the solar neighborhood. 

Subject headings: galaxies: abundances, galaxies: elliptical and lenticular, galax- 
ies: evolution, galaxies: individual: Draco dSph, galaxies: Local Group 

1. Introduction. 

Galaxies like the Draco dwarf spheroidal may be the building blocks of the Universe (Koo 
et al. 1994). Draco is among the faintest known galaxies and indeed, at 2 x lO^L©, ranks 
among the middling globular clusters in luminosity. However, its very extended distribution 
of stars {D ~ 0.5 kpc) clearly distinguishes it from the comparatively tightly condensed 
globulars. Moreover, assuming dynamical equilibrium leads one to very high estimates for 
the mass-to- light ratio (Aaronson 1983), with recent multifiber spectroscopy yielding a re- 
markable M/L = 90 ± 15 (Armandroff, Olszewski, & Pryor 1995). If the dark matter is 
primarily stellar, as might be inferred from the findings of recent microlensing experiments, 
then the Initial Mass Function (IMF) in these galaxies might have been very different from 
the present-day IMF in our own Galaxy. It is clearly important to try to piece together the 
star formation histories of such galaxies as this may have important consequences for ideas 
concerning the formation and evolution of more luminous galaxies like our own. 

The color-magnitude (CM) diagram of the Draco dwarf was first studied by Baade & 
Swopc (1961), and found to be generally similar to the CM-diagrams of metal-poor globular 
clusters with the exception of having a rather broad giant branch. Baade & Swope also 
first detected the "anomalous" Cepheids, more luminous than the period-luminosity relation 
for type II Cepheids would predict, and now known to be characteristic of dwarf galaxies. 
Another marked anomaly in the CM diagram is Draco's red horizontal branch, which is in- 
compatible with low metallicity in a Population II system (the second parameter problem). 
This phenomenon alone could be taken to suggest that Draco is ~ 2 x 10^ years younger 
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than classic, metal-poor globular clusters like M92 (Lee, Demarque, & Zinn 1994). Other 
hypotheses have been offered to explain the second parameter phenomenon (Renzini 1978), 
and an increasing body of evidence suggests that age may be no more than a minor contrib- 
utor to horizontal branch morphology (Stetson, Vandenberg, & Bolte 1996). However, the 
debate is by no means settled, and the reader is referred to Lee et al. (1994) and Chaboyer, 
Demarque, & Sarajedini (1996) for arguments in support of the age hypothesis. 

The best ground-based CCD CM diagrams of the Draco dwarf spheroidal to date have 
been presented by Stetson, Vandenberg, & McClure (1985) and Carney & Seitzer (1986, 
hereafter CS86). Like the Ursa Minor dwarf spheroidal (the 2nd-faintest dwarf spheroidal 
known), Draco appears to be as old as the oldest globular clusters, though CS86 and others 
discuss evidence from blue stragglers for an intermediate-age population as well. 

The superb resolution afforded by the Wide Field/Planetary Camera 2 (WFPC2) on 
the Hubble Space Telescope permits a remarkable improvement in our ability to study the 
stellar content of the Draco dwarf spheroidal as well as other galaxies in the Local Group. 
Here we describe the first results of an HST program to investigate the stellar population of 
this galaxy. 

2. Observations and Photometry. 

The images we presently discuss were taken on June 9, 1995 with WFPC2. The region 
selected for WFPC2 imaging is shown overlaid on the Digital Sky Survey in Figure |l] (Plate 
1), and is roughly centered on Field 1 of CS86. One 200s and two 1000s exposures were 
taken in F606W (~ V), while the F814W (~ /) images comprise 200s, 1100s, and 1300s 
integration time. F606W was chosen over F555W to give a somewhat higher signal-to-noise 
ratio. In the process of cosmic ray removal, the exposures in each filter were combined into 
single, composite images from which subsequent photometry was derived. A mosaic of the 
coadded WFPC2 frames in F814W is shown in Figure || (Plate 2). 

Photometry on the three Wide-Field (WF) detectors (0"0996 pix~^) was carried out us- 
ing a combination of aperture photometry and ALLFRAME PSF-fitting photometry (Stet- 
son 1994). ALLFRAME is the latest development in the DAOPHOT series of photometry 
packages and differs from its predecessors primarily in its ability to use information from 
many individual frames simultaneously to better constrain the final fitted magnitudes. ALL- 
FRAME is applied to the data a/i(er object detection and aperture photometry is carried out 
using DAOPHOT II (Stetson 1987). 

To push the star detection limit to as faint a level as possible, we coadded all images 
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taken through the same filter. Coadding the F606W and F814W images prior to object 
detection might in principle seem to offer an even fainter detection limit, but doing so would 
have severely compromised the reliability of subsequent completeness determinations. Two 
DAOPHOT detection passes were carried out, separated by photometry and subtraction 
of all stars found in the first pass. Faint stars hidden within the PSF skirts of brighter 
companions were thereby revealed and subsequently added to the list of detected stars to 
be measured. Point-spread functions for each detector and each filter were constructed 
within DAOPHOT using the ~ 20 brightest and most isolated point sources in each frame. 
ALLFRAME was then applied to the F606W and F814W images simultaneously. A total of 
4429 stars were ultimately detected and measured in the three WF frames. 

Due to the degree of undersampling and relatively poor PSF representation in the WF 
chips, PSF-fitting of brighter stars is subject to greater uncertainties in the absence of crowd- 
ing than is pure synthetic-aperture photometry. However, crowding and poor signal-to-noise 
ratio become ever increasing problems as one considers fainter and fainter stars, and the 
ALLFRAME-derived magnitudes yield a tighter main sequence than is obtainable using 
aperture photometry. In order to retain the small uncertainties associated with aperture 
photometry of bright, uncrowded giants and to simultaneously extend the main sequence as 
far as possible, the color-magnitude diagram was constructed using both aperture and PSF- 
fitted magnitudes. The aperture magnitudes were used whenever the predicted total light 
within an aperture of radius 2 pixels from all other resolved stars (based on the known char- 
acteristics of the WFPC2 PSF) within 2" was less than 1% of the light from the star centered 
in the aperture. Corrections to a 0"5 aperture were made in each case, and the synthetic 
zero points of Holtzman et al. (1995) were used to compute WFPC2-system magnitudes. 

Completeness tests were carried out by adding ^ 100 stars at a time to each of the 
F606W and F814W images. Stars were added in integer-F814W steps, and with F606W - 
F814W ^ 0.9 mag. The frames were then processed using DAOPHOT II and ALLFRAME 
in a manner identical to that applied to the original data. The results of these tests are 
shown in Figure ^, which shows that the 50% completeness level occurs at F814W ~ 26.2. 
The faintest stars we can reliably measure therefore have masses of ~ O.GM©. The tests also 
show that the photometric uncertainties are ±0.06 mag rms at F814W = 23, and ±0.11 mag 
rms at F814W = 25. 



3. Discussion 
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3.1. Morphology of the Color- Magnitude Diagram 

The F814W vs F606W-F814W color-magnitude (CM) diagram, constructed using both 
aperture photometry for uncrowded stars (269) and ALLFRAME fits for crowded stars 
(4156) is shown in Figure An electronic version of the photometry table is available on 
request from CJG. 

The morphology in Figure | above the turnoff is very similar to that found by Stetson 
et al. (1985) and CS86. Owing to the superior resolution afforded by the HST data, we 
are able to extend the main sequence 2 magnitudes fainter than these ground-based studies. 
The turnoff region is well resolved and reveals that the bluest main-sequence stars have 
F606W-F814W = 0.36. There are evidently fewer blue stragglers than in the CM diagram 
of CS86, consistent with a ratio of ~ 5 in field of view. Unlike CS86, we see no evidence 
for distinct multiple turnoffs; most of the stars within the WFPC2 field of view appear to 
be approximately coeval. We caution however that, while our photometry may be more 
accurate, we have sampled only one fifth as many stars in the turnoff region as have CS86, 
and it may be premature to rule out the multiple turnoffs in Draco. Additional WFPC2 
fields of Draco will be required before we can seriously address this issue. 



3.2. Age 

In Figure |^ we compare the Draco CM diagram with fiducial sequences for the Galactic 
globular clusters M68 (Walker 1994) and M92 (Heasley & Janes 1996), which have [Fe/H] 
= -2.09 and -2.24 dex, respectively (Djorgovski 1993). The Draco magnitudes have been 
dereddened assuming E(i?-V) = 0.03 (Stetson 1979a) using the absorptions tabulated by 
Holtzman et al. (1995) for stars of K5 spectral type. The globular sequences were dereddened 
and translated using the color excesses and distances tabulated by Peterson (1993). A mean 
abundance for giants in Draco has been spectroscopically measured to be [Fe/H] —1.9 dex 
(Lehnert et al. 1992). Adopting a mean horizontal branch magnitude (my) = 20.07 ± 0.03 
(Stetson 1979a), E(5 -V) = 0.03 (Stetson 1979b), and My = 0.15 [Fe/H] + 0.82 (Carney 
et al. 1992), we use (m— M)o = 19.48. The isochrones shown in Figure ^are from Vandenberg 
& Bell (1985), with colors from the V — K T^ff relations for dwarfs of Johnson (1966), 
and for giants of Ridgway et al. (1980). V — K colors have been transformed to V — I 
empirically and in the correct metallicity regime via the preliminary multimetallicity color- 
color diagrams of Worthey & Fisher (1996), and V — I have in turn been converted to 
F606W-F814W using the transformation equations in Holtzman et al. (1995). 

A sequence of normal points representing the color-magnitude sequence was generated 
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by defining an inclusion envelope around the main sequence, sub-giant, and giant branches, 
and computing mean magnitudes and colors in 0.4 magnitude- wide bands. Due to the 
scarcity of giant branch stars and the associated large uncertainties in their mean colors, 
the locus of the giant branch above F814W = 20.5 was estimated by eye. This sequence, 
resampled using spline fits, was then compared using a maximum-likelihood technique with 
isochrones ranging in age from 12 to 18 Gyrs, and in [Fe/H] ranging from -1.2 to -2.2. The 
normal points were uniformly weighted, and the best-fitting isochrone had an age of 16 Gyrs 
and [Fe/H] = -2.2. Given the weighting of the data, a formal estimate of the uncertainty 
was not straightforward, and the fitting uncertainty was instead estimated by generating one 
hundred realizations of the sequence of normal points (using Gaussian deviates based on the 
dispersion in the mean colors in each magnitude bin) and comparing them to each of the 
isochrones. This yielded an age uncertainty, based purely on random and fitting errors, of 
0.5 Gyr. 

Shown in Figure |^ are isochrones corresponding to an age of 16 Gyr and [Fe/H] ranging 
from -1.2 to -2.2 dex. An age of 16 Gyrs is not atypical of ages found for several metal-poor 
Galactic globular clusters using this set of isochrones. However, despite the small uncertainty 
due to random effects, it is apparent from Figure ^ that there are clear, systematic differences 
between the isochrones and the data. Whereas the fiducial sequences of M68 and M92 agree 
reasonably well with the CM distribution of Draco stars, the isochrones generally appear 
too red, particularly along the giant branch. This is undoubtedly due in large part to the 
roundabout technique used above to transform values of T^ff to F606W-F814W colors. 
Given the numerous calibrations involved and the systematic discrepancies which can arise, 
it is clearly dangerous to place much faith in the age we derive using this method. 

A calibration-independent method for measuring relative agesQ has been described by 
Vandenberg, Bolte, & Stetson (1990). This involves measuring the difference in color be- 
tween the turnoff and the giant branch, which is relatively insensitive to metallicity but is a 
monotonic function of age. In Figure ^ we show the fiducial sequence for Draco (generated 
as described above) compared with similar sequences for the Galactic globular clusters M68 
and M92. Also shown are isochrones having [Fe/H] = -2.0 and ages ranging from 13 to 18 
Gyr. The globular cluster sequences and the isochrones have been shifted horizontally to 
match the Draco sequence at the bluest point of the turnoff region, and vertically to agree 
with one another at a point 0.05 magnitudes redward of the turnoff on the main sequence. 
The Draco sequence evidently agrees very well with the globular cluster sequences, though 



"'^We emphasize that our conclusions refer to relative ages. If RR Lyrae stars are as luminous as claimed 
by Reid (1997) and Alcock et al. (1997) (My = +0.20), then the absolute ages of the globular clusters (and 
Draco) are reduced by ^ 3 Gyrs. 
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there are small differences in the shape of the turnoff and the slope of the giant branch. 

The relative ages between sequences are best measured in this technique by comparing 
the colors of the giant branch, where the color separation is largly independent of magni- 
tude. Looking at the isochrone spacing, it is evident that the color difference between giant 
branches is not uniform, but depends on absolute age. If we calibrate age differences using 
the oldest (leftmost) two isochrones, then the mean color difference measured in the region 
20 < F814W < 21 (below the HB and in a region where the giant branch is relatively well 
sampled) yields a relative age difference between Draco and either of the globular clusters of 
considerably less than one Gyr. Note that using increasingly younger isochrones to establish 
the color difference-age calibration would reduce the inferred age differences even more. 

Given the photometric scatter apparent on the main sequence, the presence of blue 
stragglers or large numbers of binary stars could conceivably affect the determination of the 
mean colors of the turnoff and the main sequence. To estimate our sensitivity to these effects, 
we have recomputed the Draco normal points using an inclusion envelope alternately more 
exclusive of stars on the blue side of the turnoff and on the red side of the main sequence. 
In the latter case, we have excluded approximately three times as many outliers to the red 
of the main sequence as to the blue. Using these computed normal sequences (in addition 
to a number of hand-sketched realizations) and bringing all other sequences into alignment 
as before, we estimate our measuring uncertainty to be ±1 Gyr. The distribution of relative 
ages measured in this way is slightly asymmetric in the sense that an older Draco is preferred 
over a younger. 

Additional uncertainty arises from possible errors in the color terms used to transform 
between photometric systems. Two such color terms are used here: the term used to trans- 
form the ground-based data into V-I, and that used to transform V-I to F606W-F814W. The 
close similarity in age we infer for M68 and M92 agrees with the results of Vandenberg et al. 
(1990) based on independent B, V data sets and suggests that errors in the ground-based 
color terms are likely to be small. Errors in the transformation from V-I to F606W-F814W 
arise from inaccuracies in the fitting formula used on the synthetic photometry results and 
from the error made by applying a transformation based on stars of solar metallicity to a 
metal-poor population. We estimate that fitting errors would introduce an error of order 
0.02 magnitudes in the difference in color between the turnoff and the giant branch. This 
translates to an age uncertainty of about 2.3 Gyrs. 

Applying transformation coefficients based on solar-metallicity stars to the more metal- 
poor stars of M68 and M92 has the effect of making the globular cluster giant branches in 
Figure |^ appear somewhat bluer relative to Draco than they should. Using Kurucz model 
atmospheres, we find that stars at the turnoff and on the giant branch are fainter in F606W 
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by 0.014 and 0.024 magnitudes, respectively, than what we compute using the Holtzman 
et al. (1995) coefficients. The F814W transformation for these stars is much less sensitive to 
metallicity, and we estimate that the differences in F606W-F814W color between the turnoffs 
and the giant branches of M68 and M92 should be about 0.014 magnitudes larger than shown. 
Correcting our measured result accordingly and combining the two major sources of error 
above, we conclude that Draco is 1.6 ± 2.5 Gyrs older than either M68 or M92. 

Given the very different horizontal branch morphologies of Draco and the two globular 
clusters, what can we say about the 2nd parameter problem? In Figure ^ we partially repro- 
duce the [Fe/H] vs. HB morphology diagram of Lee et al. (1994). We have determined the 
morphological type quantity (B-R) / (B+V+R) for Draco using the color-magnitude diagrams 
of both Stetson (1979a) and CS86. Counting stars on the blue and red sides of the insta- 
bility strip as well as likely RR Lyrae variables, we find (B-R)/(B+V+R) = -.77 and -0.70 
for the two data sets, respectively. We plot the mean of these two values in Figure ^ along 
with similar quantities tabulated by Lee et al. (1994) for M68 and M92. For both uniform 
mass-loss models and models in which the mass-loss rate increases with age (Reimers 1977), 
the Draco point appears to fall about 2 Gyr below the isochrone which best matches M68 
and M92. As this is only 1.4x our estimated age uncertainty, we cannot completely rule 
out consistency between Draco and an age-dependent HB morphology of the type modeled 
by Lee et al. . WFPC2 photometry of old globular clusters, sufficient to reduce the relative 
age uncertainty to ~ 0.5 Gyr, will be required before we can put ffimer constraints on the 
models. 



4. Star Formation History 

Zinn (1978) and Lehnert et al. (1992) have demonstrated an abundance spread of ap- 
proximately 1 dex in Draco. Aside from this chemical inhomogeneity, how well does the 
Draco dwarf spheroidal galaxy fit the simple stellar population paradigm of Renzini & Buz- 
zoni (1986)? We have seen that the turnoff CM diagram morphology of Draco is very similar 
to that of M92; globular clusters are the prototype simple stellar populations, coeval sys- 
tems of initially homogeneous composition. Clearly, Draco is predominantly a simple stellar 
population originating as early as anything in the Galaxy. 

Close examination of the CM diagram shows a significant population of blue stragglers 
(see also CS86). Such stars are often seen in simple stellar populations, including both an- 
cient globular and younger open star clusters, but they can also result from an extended 
star formation history. We can count them in color profiles (Figure formed out of half- 
magnitude strips from the turnoff region. The surface density of foreground stars in this 
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region of the CM diagram (specifically a region 0.3 magnitudes wide in color and 2 magni- 
tudes high in F606W) is predicted by the Galactic model of Bahcall and Soneira (1980) to 
be 17 per square degree. Thus, there is only a 2% probabihty of finding a single foreground 
star among the ~ 17 stars we identify as blue stragglers (roughly speaking, those stars with 
F814W < 23 and F606W-F814W < 0.3). There are between 60 and 70 stars in the same 
region (V < 23 and B — V < 0.3) of the CM diagram of CS86. Given that the area surveyed 
by CS86 is about 5 times that of the area covered by our WFPC2 field, a ratio of 3.8 ± 1 in 
the number of blue stragglers indicates reasonable agreement. 

Simulations of the turnoff region using the Yale isochrones (Green, Demarque &: King 
1984) require 6% of the mass to be contained in a second population, half the age of the 
primary population, in order to fully account for the blue stragglers in Figure ^. We note 
that 6% is an upper limit on the size of a second younger component of Draco, since blue 
stragglers can be an intrinsic part of a simple stellar population. 

A putative extended star formation history strengthens the link between Draco and more 
fully developed examples of this phenomenon, such as Carina (Smecker-Hane et al. 1994; 
1996). The difference between Draco and Carina and galaxies such as LGS 3 (Mould 1997) 
which are still forming stars today (at 10% of the original rate in the latter case) is then simply 
one of degree. Draco simply had a shorter tail to its initial burst of star formation. There is 
considerable evidence that wind-induced mass-loss cannot account for the stellar populations 
seen in several Local Group dwarfs, and that stripping or capture of gas must have occurred 
(see Gallagher & Wyse 1994 for a review). Draco is currently among the nearest of the dwarf 
galaxies and, if its orbit is largely contained within its present distance, the probability of 
capturing significant amounts of rejuvenating circumgalactic gas may be substantially smaller 
than for dwarfs in the outskirts of the Local Group. Alternatively, a high ultraviolet flux 
from the precursors to the Galaxy and M31 may have removed substantially more gas from 
the nearby dwarfs (van den Bergh 1994). 

None of these remarks bear at all on the better-known curious phenomenon of Draco's 
red horizontal branch. The finding that Draco is as old or older than M92 requires us to 
focus on the problem of understanding Draco's RHB. Higher S/N (and preferably F555W- 
F814W) photometry performed differentially with respect to [Fe/H] = -2 clusters have the 
potential to sharpen our conclusions further. If the relative age uncertainty could be reduced 
to 0.5 Gyr, one would have to turn to parameters such as cluster density (e.g. Fusi Pecci 
et al. 1996) as the second parameter in the globular cluster system. The age spread of the 
halo, which is the defining feature of the Searle & Zinn (1978) galaxy formation scenario (as 
opposed to Eggen, Lynden-Bell, & Sandage 1962) would then need to be defined by the ages 
of the low density element of the halo, namely Draco, its sibling dwarf spheroidals, and the 
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Galactic bulge. 

4.1. The Initial Mass Function 

Figure |^ shows the differential luminosity function (LF) for the stars in our sample. Also 
shown is a LF computed from Walker's (1994) photometry of M68. No completeness correc- 
tions have been applied to the M68 data, and the counts have been arbitrarily normalized 
to our data at F814W = 23. Note that the number of measured stars in the M68 sample is 
about 6x larger than in the Draco sample. The Draco LF appears to be somewhat steeper 
on the giant branch, and somewhat shallower below the turnoff. Uncertainties introduced 
by differences in noise characteristics, resolution, and detection efficiency make it difficult to 
draw firm conclusions but, over the (admittedly narrow) mass range considered here, there 
does not appear to be an obvious bias towards lower-mass stars which might be held to 
account for the very high estimates of Draco's mass-to-light ratio. Based on the present 
data set, all we can say is that if the bulk of the matter is in the form of stars, such stars 
must have masses less than O.GM©. 

The observed LF can be used to place limits on the initial mass function (IMF) in Draco. 
The current observations are of stars with masses between 0.6 and 0.9 Mq. We parameterize 
the IMF in this mass range by a simple power law, with dN/dM oc M°. Even stars of this 
low mass evolve significantly over a Hubble time, so the inferred IMF slope depends on the 
age of the stellar population in Draco. 

We construct model LFs using a variety of IMFs and compare each model distribution 
with the observed one to see if it is plausible that they are drawn from the same parent 
population. The construction of model LFs allows us to incorporate photometric errors 
(both random and systematic) and incompleteness into the analysis; the incompleteness 
and errors which we apply to the model LFs are those determined from our artificial star 
experiments. We can also include the effect of a binary star contribution if we make an 
assumption about the frequency and distribution of the relative masses of stars in binary 
pairs. For Draco, we assume a metallicity of [Fe/H]=-2, a distance and extinction as quoted 
above, and use Vandenberg (1983) isochrones (from which Worthey's models are derived for 
low mass stars) along with Kurucz model atmospheres to convert from masses to assumed 
luminosities. 

We define as unacceptable IMF slopes for which the model LF and the observed LF are 
inconsistent with each other at the 95% confidence level as determined by a Kolmogorov- 
Smirnoff test. With this definition, we constrain the IMF slope to fall within the range 
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a > —1.00 for an assumed age of 18 Gyr, 
— 1.52 > a > —1.98 for an assumed age of 15 Gyr, 
-2.14 > a> -2.30 for an age of 12 Gyr, and 
-2.48 >a> -2.56 for an age of 10 Gyr 

when we fit the F606W LF over the range 1 < Mpeoew < 7.25 with stellar models of a 
single metaUicity, [Fe/H]=-2. The dependence on age is easy to understand: as a population 
evolves, the LF around the turnoff steepens, leading to a flatter inferred IMF for an older 
population. These results were derived assuming a binary fraction of 0.5 with uncorrelated 
masses of the binary stars; with no binaries, we get slopes which are slightly shallower. 
Similarly, small variations on the above allowed ranges are derived with different assumptions 
about the mctallicity, distance, and/or a different adopted magnitude range over which the 
comparison with the observed LF is made. 

If we fit the F814W LF, we infer a range of possible IMF slopes which are slightly 
steeper than those inferred from the F606W LF. For no binaries, the inferred slope range 
overlaps between the two bandpasses; for a binary fraction of 0.5, all of the acceptable IMF 
slopes for F814W are steeper than those allowed by the F606W LF. 

For comparison, the IMF in this mass range in the solar neighborhood is proportional 
to M~^'^ (Kroupa, Tout, & Gilmore 1993). The allowed slopes for Draco are remarkably 
similar to this value. 

5. Conclusions 

Based on photometry of WFPC2 images of the Draco dwarf spheroidal, we conclude the 
following: 

• Draco is predominantly an old system, being 1.6 ± 2.5 Gyrs older than either M68 or 
M92. 

• We see no evidence for multiple turnoffs which would be indicative of episodic star 
formation. 

• The presence of significant numbers of blue stragglers may indicate an extended period 
of star formation, though the numbers of younger stars < 6% of the first generation. 
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• The luminosity function looks in all respects nearly identical to that of the old globular 
cluster M68. 

• The slope of the initial mass function appears to be very similar to that of the IMF 
determined for the solar neighborhood. 

Future work should be aimed at refining the relative ages of Draco and blue horizontal 
branch clusters hke M92. Obtaining both higher signal-to-noise ratio data as well as WFPC 
color-magnitude diagrams of metal poor globular clusters should allow us to reduce the 
relative age uncertainty to ~ 0.5 Gyr. Such data would considerably strengthen constraints 
on the process and mechanisms which led to the formation of our Galaxy. 

We are grateful to Gary da Costa for a critical reading of an earlier version of this 
manuscript. This research was conducted by the WFPC2 Investigation Definition Team, 
supported in part by NASA Grant No. NAS7-1260. 
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Fig. 1. — Digitized Sky Survey image of the field containing the Draco dwarf spheroidal. The 
outlines indicate the HST WFPC2 field and Fields 1 and 2 of Carney and Seitzer (1986). 
The entire field shown subtends 15' on a side. 

Fig. 2. — Mosaic of the coadded F814W frames. 

Fig. 3. — Completeness fraction as a function of F606W and F814W magnitude. Complete- 
ness tests were made using artificial stars of color F606W - F814W a; 0.9 mag. 

Fig. 4. — HST color-magnitude diagram of the Draco dwarf spheroidal. Most of the stars 
on the giant branch were measured using aperture photometry, while the majority of stars 
below the turnoff were measured by fitting PSFs. 

Fig. 5. — The color- magnitude diagram of the Draco dwarf spheroidal compared with fiducial 
sequences for the metal-poor globular clusters M68 and M92. Superimposed are Worthey 
models of age 16 Gyrs and spanning metaUicities of [Fe/H] = -2.2 on the left to [Fe/H] = 
-1.2 on the right. The globular cluster seqiicnccs have been transformed from V and / to 
F606W and F814W using the transformation coefficients given by Holtzman et al. 1995. 

Fig. 6. — The fiducial sequence for Draco compared with similar sequences for the Galactic 
globular clusters M68 and M92, using the relative-age-measuring technique described by 
Vandenberg, Bolte, & Stetson (1990). The M68 and M92 sequences, measured in V and /, 
have been transformed to F606W, F814W using the synthetic transformation discussed in 
the text. The cluster sequences and the Worthey models (having [Fe/H] = -2.0 and ages 
ranging from 13 Gyr on the right to 18 Gyr on the left) have been shifted horizontally to 
match the Draco sequence at the bluest point of the turnoff region, and vertically to agree 
with one another at a point 0.05 magnitudes redward of the turnoff on the main sequence. 
The filled circles correspond to the locus of the giant branch determined as described in the 
text. 

Fig. 7. — A partial reproduction of Lee et al. 's (1994) Figure 9a, which they use to demon- 
strate a relation between horizontal branch morphology and age. The solid lines correspond 
to models which assume uniform stellar mass loss and are separated by 2 Gyrs. The dashed 
fines are also separated by 2 Gyrs and assume a mass-loss rate which increases with age. The 
plotted values for various globular clusters are taken from the compilation of Lee et al. (1994), 
and the HB type for Draco was computed using the color-magnitude diagrams of Stetson 
(1979a) and Carney & Seitzer (1986). 

Fig. 8. — Color profiles generated by counting stars lying within half-magnitude-wide F814W 
bins. We attribute the excess number of blue stars over a putative Gaussian color distribution 
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to the presence of ~ 17 blue stragglers. 

Fig. 9. — The differential luminosity function of the stars in Draco. The error bars reflect 
Poisson errors only, and the dashed histogram shows the effect of dividing by the complete- 
ness function. Note that no account has been taken of bin-jumping due to photometric 
errors. The dotted hne shows a luminosity function derived from Walker's M68 photometry 
without having applied completeness corrections. 



This figure "figl.gif" is available in "gif" format from: 
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http://arXiv.org/ps/astro-ph/9709259vl 




.2 .4 .6 .8 1 1.2 

F606W- F814W 




.4 .6 .8 

F606W- F814W 



20 
15 
10 
5 

20 
15 
10 




40 
30 
20 
10 



22.0 



22.5 



Q 1_J I \ \ \ I \ \ \ \ \ I I I \ \ \ I L 

.2 .4 .6 .8 1 

F606W- F814W 



